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1. Introduction 



1.1. Background 

The theory of electronic transport in macroscopic normal metals has traditionally 
employed the Boltzmann or quasiclassical approach to determine an electronic 
distribution function in the presence of impurity and phonon scattering, electric 
and magnetic fields, and band-structure effects [ 1 ]. From this distribution function 
the spatially-averaged current is calculated and conductivity. Hall conductivity 
and non-linear I-V characteristics are obtained. Successive scatterings are treated 
as independent and results are averaged over configurations of impurities. Sample 
geometry is normally not considered relevant as the macroscopic conductivity is 
an intensive quantity. 

Beginning with the prediction and discovery of the weak localization effect 
in metals [2-4] it was appreciated that at sufficiently low temperatures even in 
macroscopic conductors the assumption of independent scatterings misses qual- 
itatively new effects related to electronic phase-coherence. However with the 
discovery of the hc/e period Aharonov-Bohm effect in single metal rings [5,6] 
and reproducible mesoscopic conductance fluctuations [7-14] it became clear that 
transport in conductors on the size scale of a micron or below would exhibit prop- 
erties very different from those of macroscopic conductors. As the physics of 
these systems became better understood it was recognized that these properties 
are not characteristic of a fixed size scale but rather depend on the temperat- 
ure, T, fermi wavelength, kf and degree of disorder in the material. None of 
the phenomena discussed in this volume occur in micron-size conductors at room 
temperature; conversely mesoscopic phenomena such as conductance fluctuations 
have been observed in conductors of size up to 100 microns at milliKelvin tem- 
peratures in samples of high purity and small kf. Hence the term mesoscopic was 
introduced to denote a physical regime for transport and not a specific size scale. 

1.2. Novel properties of mesoscopic conductors 

The novel transport properties of mesoscopic conductors may be roughly sum- 
marized as follows: 
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Fig. 1 . Measurement of non-local Aharonov-Bohm oscillations in conductance by Umbach et al. [17]. 
For sample with loop attached outside the current path the periodic A-B oscillations are superimposed 
on the aperiodic universal conductance fluctuations. Adapted from [18]. 

Sample-specificity. Mesoscopic conductors of the same material fabricated in 
the same manner will have different measured transport coefficients. This is il- 
lustrated by the different (but reproducible) conductance fluctuations as a function 
of magnetic field (or other external parameters) measured in mesoscopic samples 
(see Fig. 1). 

Non-locality. The measured transport properties of a mesoscopic conductor do 
not depend solely on the portion of the sample between the measuring points. If 
the phase-coherence length [15,16] is sufficiently large then regions of the sample 
outside the current path strongly influence the measured properties. 

This is illustrated dramatically by the experimental data shown in Fig. 1 [17]. 
In this experiment a mesoscopic ring was attached to a conducting microstructure 
outside the current path and was found to a yield a measurable hc/e Aharonov- 
Bohm oscillation in the magnetoresistance due to the interference between elec- 
tronic trajectories along the main current path and those which detour around the 
ring. 
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Fig. 2. Magnetoconductance measurements taken for three different configurations of current and 
voltage leads. The AB oscillations have been filtered out of the data. Note the complete lack of 
symmetry around B = for a given measurement, but the mirror symmetry of the top and bottom 
traces. This symmetry is required by reciprocity relations to be discussed below. Reprinted from [6]. 



Violation of macroscopic symmetries. A number of symmetries found in the 
properties of macroscopic conductors are violated by mesoscopic conductors [11]. 
These are symmetry properties that do not follow immediately from microscopic 
symmetries but also require further conditions (e.g. symmetry of the scattering 
potentials on average) which were often not explicitly recognized until they were 
found to be violated in mesoscopic conductors. One striking example is the lon- 
gitudinal magnetoresistance. A standard measurement of R xx involves injecting 
and collecting a current at two ends of the sample and measuring the voltage 
induced between additional probes attached along the sample with the second 
voltage probe displaced a distance much larger than the width of the wire. Under 
these conditions in a homogeneous macroscopic conductor the Hall contribution 
to the measured resistance is negligible and one measures a magnetoresistance 
which is symmetric in B. However this symmetry does not follow from the funda- 
mental Onsager symmetries of the non-local resistivity tensor (see below); it also 
relies on approximate symmetries of the transmission probabilities for electrons 
which are violated in the mesoscopic regime. In fact it has been shown [20,19,21] 
that no matter how one attaches the voltage probes in a phase-coherent sample 
the fluctuating part of the measured resistance arises roughly equally from the 
symmetric ("longitudinal") and anti-symmetric ("Hall") resistivities so that the 
resulting magnetoresistance has no symmetry under field reversal (see Fig. 2). 
Similarly the I-V characteristic in mesoscopic samples becomes asymmetric at 
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Fig. 3. Weak localization effect in 3 cases: Top trace is a wire much longer than the phase coherence 
length. Middle and bottom traces are from similar mesoscopic wires attached to narrow and wide 
probes respectively. Reprinted from [23]. 



extremely small voltages corresponding to internal electric fields much smaller 
than are needed to observe non-linearities in conventional conductors [22]. 

Dependence on measurement geometry. Measured resistances also depend on 
what type of probes are employed. In Fig. 3 is shown a comparison [23] of 
weak localization magnetoresistance for two mesoscopic samples which differ 
primarily in the nature of the voltage and current probes used. The sample with 
narrow probes shows a WL amplitude roughly one order of magnitude larger than 
the comparable sample with wide probes. Intuitively this is related to the difficulty 
an electron has in returning phase-coherently to the sample from the wide leads. 
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Fig. 4. Hall resistance of a GaAs quantum wire junction with flattened corners (inset) which leads 
to inverted Hall resistance (solid line); data from a nominally square-cornered junction (dashed line) 
show quenching (vanishing slope near B = 0) but not inversion. Inset shows classical rebound 
mechanism for quenching and inversion discussed in section 4.4. Adapted from [24]. 



Dependence on sample geometry. Much of this chapter will focus on extremely 
high-mobility semiconductor microstructures (as opposed to disordered metallic 
conductors which are treated in detail in chapter [15,16]). In the high-mobility 
microstructures the dominant source of resistance is the sample geometry itself. 
This is illustrated in Fig. 4 where we show Hall resistance data [24] from GaAs 
quantum wire junctions with differing geometries. 

These results taken together demonstrate that the concept of resistivity as an 
intensive quantity breaks down at the mesoscopic scale. Instead one needs to de- 
velop a transport theory based on sample-specific quantities: i.e. quantities which 
are meaningful for a given sample measured in a given manner. We now review 
the main approaches which have been developed to describe mesoscopic transport. 
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2. Sample-specific theory 

2.1. Conductance coefficients 

A typical mesoscopic transport measurement is performed by imposing a very 
low frequency AC current through the microstructure and measuring the voltages 
induced between various auxiliary probes (this is called a multiprobe resistance 
measurement). At sufficiently small currents there must be a linear relationship 
of the form 



where I m is the current measured in probe m and V n — Vq is the difference between 
the voltage measured at probe n and a reference voltage V . The g mn are the 
conductance coefficients and are independent of voltage at small currents. In a 
typical measurement a current I is injected in a given probe (probe 1) and ex- 
tracted through probe 2 and the other "voltage probes" draw no current; hence 
Eq. (2.1) is to be solved with I\ = /, I2 = — /, I m — 0, m — 3,4, . . .. The 
resistances are given by R\2,mn = (V m — V n )/L where the first pair of in- 
dices denote the current leads and the second the voltage leads. Since the res- 
istances must be independent of the reference voltage Vo and the currents must 
satisfy Kirchhoff's law (X^m=i = 0) the conductance matrix g must satisfy 
]Cm=i 9 mn = X^n=i 9mn = 0. These two conditions are independent in general 
as g need not be symmetric (we shall find it is symmetric only in the absence of a 
magnetic field). These conditions imply that g is only invertible in the subspace 
perpendicular to the vector (1,1,1,...) and solutions of Eq. (2.1) have a non- 
uniqueness corresponding to the redefinition of zero voltage. The conductance 
matrix g contains in principle all the sample-specific information relating to the 
linear resistance of mesoscopic conductors. The challenge then is to calculate 
the g mn in a manner which includes the effects of phase-coherence and sample 
and measuring geometry which we have seen above are important. The main ap- 
proach to meeting this challenge has been to relate the conductance coefficients 
to an appropriate quantum-mechanical scattering problem. This was achieved by 
arguments due to Landauer, Imry and Biittiker; it has also now been understood 
that the same conclusion can be reached by applying quantum linear response 
theory. 

2.2. The scattering approach 




(2.1) 



71=1 



I will not provide a full historical review of the literature on this topic here, as I 
have done so elsewhere [25]. Briefly, long before the birth of mesoscopic phys- 
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ics Rolf Landauer [26] gave a physical argument leading to the conclusion that 
a truly one-dimensional phase-coherent quantum-mechanical conductor should 
have a conductance given by g = (e 2 /h)T/R, where T, R referred to the total 
transmission and reflection coefficients of the conductor(I will neglect spin except 
when needed for comparison with experiment). This formula had the intuitively 
appealing property that as the sample transparency approached unity (T — > 1) 
its resistance vanished. The first attempts to derive this formula from linear re- 
sponse theory however yielded the result g — (e 2 /h)T [27] or its many-channel 
generalization [28] g = (e 2 /h) X^t=i T a b, where a, 6 index the N c channels 
injecting from the leads. Engquist and Anderson [29] showed quite early that the 
Landauer formula really was meant to describe a four- probe measurement and 
Imry [30] subsequently understood that g = (e 2 /h)T described an ideal two- 
probe measurement. He pointed out that the non-vanishing resistance predicted 
by this formula for a "perfect conductor" (T — ► 1) actually represents an ideal 
quantum contact resistance which is unavoidable in a two-probe configuration. 
Very shortly thereafter this quantized contact resistance was directly observed in 
semiconductor point contacts in GaAs heterostructures [31,32]. Biittiker [33,34] 
completed the development of the scattering approach by focusing on the en- 
tire conductance matrix and showing that a generalization of Landauer's original 
"counting argument" leads to a simple and useful theory for the quantum conduct- 
ance matrix with the correct time-reversal symmetry properties. I now reproduce 
this argument. 

The sample whose resistances are being measured is assumed to be attached 
to Ni perfect phase-randomizing reservoirs all at chemical potential /n through 
ideal leads which can serve either as current or voltage probes. When voltages 
or currents are imposed between the reservoirs (which roughly correspond to 
massive contacts in real experiments) in general current will flow transiently into 
or out of all the reservoirs. However the chemical potentials of those reservoirs 
chosen to be voltage probes will be adjusted until these currents are cancelled, 
leaving a steady-state voltage V n on all the probes. Denote the current into probe 
m due to the voltage (measured from equilibrium) on probe n as I mn where I m = 
En=i Ann- Assuming m ^ n we can write I mn = Y^b=\ Ea=i T mn.ab- 
where I® b is the current injected in channel b from probe n in the energy inter- 
val eV n — fj.Q and T mn ,«6 is the transmission coefficient from lead n, channel b to 
lead 77i, channel a. Since the reservoirs are assumed to feed each incident channel 
equally, the current I® b = evb{dn b jde)eV n where Vb is the velocity in channel b 
and dn b /de is the density of ingoing states in channel b in lead n (we assume all 
leads have the same channel states for simplicity; this assumption is not required). 
If one imagines the lead to be perfectly translationally invariant in the longitud- 
inal direction then one has the equality Vbdn* /de — 1/2-Kh — l/k. This identity 
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holds for non-interacting quantum particles with arbitrary transverse confinement 
in the presence of a uniform magnetic field of arbitrary strength. The transla- 
tional invariance in the longitudinal direction means that the propagating states 
may be labelled by a wavevector k and a mode index b such that their longitudinal 
velocity is hv b = de b (k)/dk = (27rdn^"/de) _1 . Because this only holds for non- 
interacting particles the scattering approach only holds for quasi-particles which 
cany unit charge as in fermi liquids and must be non-trivially generalized to treat 
mesoscopic superconductors or the fractional quantum Hall regime. Combining 
these results yields 

N L 2 N L N l 

Ira = Yl I mn = ^ ^(^n " N c 6 mn )V n = £ mn V» (2.2) 

n=l n=l n=l 

(the term N c 6 mn arises when the argument is applied to the case m = n and we 
have set V = 0). Here 

Tmn — ^ ^ Tmn.ab Tl {t mn tj nn } _ , (2.3) 
a,b=l 

where t mn is the transmission amplitude matrix between the channel states and 
tmm is to be interpreted at the reflection matrix. One can easily see that the 
unitarity relations for the S-matrix ensure that Eq. (2.2) is consistent with the 
conditions on the conductance matrix Sm=i 9 mn = 12n=i 9 mn = and do not 
alone require that g be symmetric. 

2.3. Time-reversal symmetry constraints on g 

Since Eq. (2.2) directly relates the g mn to a current-conserving quantum S-matrix, 
it is straightforward to find the constraints time-reversal (TR) symmetry imposes 
on the g mn (B) and hence the magnetoresistances. S relates the two N c Nl- 
dimensional vectors of incident and outgoing fluxes, I, O: S(B)I — O. Neg- 
lecting spin, TR symmetry implies that if the magnetic field is reversed then initial 
and final states may be interchanged (as long as they are complex conjugated): 
S(-B)0* — I*. Using the unitarity of the S-matrix this equation implies 

S T (-B) = S(B). (2.4) 

When B — this immediately implies that g is a symmetric matrix. An important 
case is if the conductor is only attached to two-probes, so that the conductance 
matrix is only 2x2. Then the requirements that its rows and columns add up to 
zero imply that the matrix is always symmetric, independent of the presence of a 
magnetic field; hence that gu = 521 = 9 is symmetric in magnetic field. In fact 
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two-probe measurements of mesoscopic conductance fluctuations do find that the 
complicated sample-specific magnetofingerprint is symmetric in B. 

In the general multiprobe case the conductance formula (2.2) and time-reversal 
symmetry of the S-matrix only imply the following symmetry: 

g mn {B) = g nm {-B). (2.5) 

This relation does not imply any symmetry between a single resistance measure- 
ment at positive and negative fields. The symmetry under field reversal observed 
in standard longitudinal magnetoresistance measurements (and antisymmetry in 
Hall measurements) arises because the impurity potentials involved are macro- 
scopically homogeneous (the violation of these symmetries in macroscopically 
inhomogeneous conductors was known previously [35]). Although these sym- 
metries are often referred to as Onsager relations we shall see that it is only the 
Onsager relation for the non-local conductivity (and the reciprocity symmetry to 
be discussed immediately below) which follows from TR symmetry alone. In 
mesoscopic transport, where quantum interference effects are known to be ex- 
tremely sensitive to the details of the impurity potential or sample geometry, typ- 
ically the spatial potentials or the measuring configuration has some asymmetry 
and one observes significantly asymmetric magnetoconducance as seen in Fig. 2. 
However Eq. (2.5) does imply an experimentally verifiable reciprocity relation 
[34]: 

^mn,m'n'(5) = Rm' n' ,mn{ — B). (2-6) 

Hence the same resistance is obtained from two different measurements: in the 
first measurement current is run from lead m to n and voltage measured between 
777,', n' at field B; in the second measurement current and voltage leads are inter- 
changed and the field is reversed. The validity of Eq. (2.6) can be seen by the 
following argument. As noted in section 2.1 above, although g is not invertible in 
general, it is invertible in the physical subspace of vectors satisfying Kirchhoff 's 
Laws. Denote the inverse r, and define unit vectors e n , n = 1,2,... Nl corres- 
ponding to unit current (or voltage) on each probe. It follows then from I = gV 
that the resistance R m n,m>n'{B) = (e m < -e n >)-r- (e m - e„). R m 'n' ,mn{-B) is 
the transpose of this with r(B) — > r(-B). But this is just equal to R m n,m'n'(B) 
as long as r satisfies r T {-D) = r(B) which is straightforward to confirm from 
Eq. (2.5). 

2.4. Resistance quantization and negative resistance 

As noted above, the general properties of the conductance matrix g and hence 
of the resistances i? m ' n ',mn follow from Kirchoff's laws and time-reversal sym- 
metry. The reciprocity symmetry Eq. (2.6) was known previously and derived 
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under classical assumptions of a local conductivity obeying Onsager's relations 
[35]. (Buttiker's argument shows that they hold even in the mesoscopic regime 
when a local conductivity cannot be defined.) The basic novelty of the scat- 
tering approach (as compared e.g. with quasiclassical transport) is the explicit 
appearance of the quantum unit of conductance e 2 /h in the conductance matrix. 
This implies that under conditions in which the transmission probabilities become 
quantized, the resistances will measure this quantum unit of conductance. Two 
experimentally realized examples are semiconductor point contacts and quantum 
Hall resistors. The semiconductor point contacts were created by dividing a GaAs 
2DEG with a split gate [31,32]. When negatively charged this gate created a con- 
striction with a width dependent on the magnitude of the gate voltage and a length 
comparable to the gate width of roughly 2500 A. In these high-mobility structures 
the electrons were negligibly back-scattered in traversing the constriction. Thus 
the structure was similar to the "perfect" conductor connecting two reservoirs 
(represented by the wide regions of the 2DEG). The transmission coefficients for 
the propagating channels in the constriction were essentially unity; the number 
of propagating channels depended on the voltage on the constricting gate lead- 
ing to conductance jumps of 2e 2 /h with decreasing (negative) gate voltage every 
time a new channel began to propagate (see Fig. 5, the factor of two is for spin). 




-2 -l.B -1.6 -1.4 -1.2 -1 



GATE VOLTAGE (V) 

Fig. 5. Conductance of a semiconductor point contact as a function of gate voltage which controls 
the size of the constriction; inset shows schematic of the structure. Adapted from [31]. 
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Although this effect follows straightforwardly from the two-probe conductance 
formula g — {e 2 /h) Ylal^b if one assumes an ideal connection between the 
"reservoir" and the "sample", its clear appearance in this experiment was not an- 
ticipated since no special effort was made to match adiabatically [36] the wide 
and narrow regions. Later theoretical work [37] showed that even in the limit 
of a completely abrupt interface which generates strong inter-channel scattering 
rather sharp transmission steps occur at T — 0. 

A second example is the integer quantum Hall effect. This is not a mesoscopic 
effect and does not require samples small with respect to either the elastic mean 
free path or the inelastic mean free path. Nonetheless it has been possible to 
develop a scattering approach to the IQHE in which the channels refer to edge- 
states at the fermi energy (see section 3.3 below) which yields some insight into 
the origin of the quantization and some predictions regarding experiments in in- 
homogeneous geometries [38,40,39]. Space limitations do not allow me to review 
this literature here. 

A final point is that the four-probe resistances obtained from Eq. (2.5) are not 
a measure of the power dissipated in the sample and are not necessarily positive. 
The power expended in transfering elctrons from the source to sink is equal to 
the product IV, but need not be dissipated anywhere in the sample, and in the 
absence of significant inelastic scattering will all appear in the sink. In contrast to 
our intuition from macroscopic conductors, any particular four-probe resistance 
can turn out to be positive or negative (with an absolute value bounded by the 
two-probe resistance) depending on the details of the particular scattering config- 
uration within the sample. Such negative resistances have actually been measured 
in microstructures, arising either due to quantum interference effects [41] or due 
to bends in ballistic microstructures [42,43]. 

3. Sample-specific linear response theory 

3.1. Model hamiltonian for mesoscopic conductors 

Equation (2.2), which relates the conductance coefficients of a mesoscopic system 
to a quantum scattering problem, has proved very useful in the theoretical descrip- 
tion of experiments on both disordered and ballistic microstructures. However 
its derivation involves the concept of an ideal electron reservoir and depends on 
the assumption that the current is carried by non-interacting fermions [44]. It has 
therefore been of interest to connect [25] this more heuristic treatment of quantum 
conductance to conventional quantum linear response theory (LRT) pioneered by 
Kubo and Greenwood [45]. The latter approach is used in most standard texts (for 
example, [46]) and is the basic technique for introducing many-body effects into 
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the theory of the macroscopic transport coefficients. Moreover any actual micro- 
scopic calculation using the scattering approach needs to introduce a hamiltonian 
since transmission coefficients are only defined for an infinite hamiltonian system 
with an asymptotic region in which scattering is absent. We will now show that 
there is a well-defined hamiltonian system for which one can derive Eq. (2.2) 
for the conductance coefficients. From the technical point of view the most in- 
teresting aspect of the derivation is that it is valid in arbitrary magnetic field and 
hence provides a rigorous formulation for calculating microscopically the Hall 
resistance (to all orders in the magnetic field). Until Buttiker's work [33] no for- 
mulation of comparable simplicity was known, and indeed there was confusion 
in the literature as to whether the Hall resistance could be expressed as a fermi 
surface property. The derivation presented here is based on that of references 



Our model consists of a finite region (the "sample") in which there is a uni- 
form magnetic field of arbitrary strength, an arbitrary static potential U(r) and 
(initially) an arbitrary electron-electron interaction V(r — r'), connected to Nl 
infinite "leads" which emanate radially outward from the sample and confine the 
electrons to non-overlapping spatial channels in the asymptotic region. The leads 
are "perfect" in the sense that they are translationally invariant in the longitudinal 
direction and the impurity and interaction potentials are assumed to vanish in the 
leads. The magnetic field is not assumed to vanish in the leads. Since the system 
is infinite, solutions to the Schrodinger equation exist at all energies above the 
threshold for propagation at infinity (which we will take to be the zero of energy). 
Since the system is described by a hermitian hamiltonian, current and energy are 
exactly conserved and in the equilibrium state no net current flows through the 
sample (although for B ^ there will be a circulating current giving rise to orbital 
magnetization). Now we wish to introduce a perturbation which represents the ef- 
fect of imposing a voltage across the sample. In a real conductor placed between 
the terminals of a battery, charge carriers are literally injected and collected at the 
terminals by means of a chemical reaction which causes charge transfer against 
the local electric field [48]. It is not obvious how to represent such a process in 
terms of a simple perturbing operator acting on the electrons. However we attempt 
to mimic the effect of the battery by imposing an external potential 0(r) cos cute 61 
where (p(r) tends to a fixed set of constants {V n }, n = 1,2,..., Nl on the leads 
(and t < 0). The external potential must be turned on adiabatically in order to 
connect it with the equilibrium state in the past [46], hence the perturbation takes 
the form: 



[40,47]. 




(3.1) 



where p is the electron density operator. 
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Note here that <f)(r) describes the external potential and its gradient gives only 
the external electric field; in the presence of this perturbation there will be an 
induced charge to linear order in the external potential which (when we include 
electron-electron interactions) will yield an additional induced potential. Hence 
the total electric field in the sample differs from — V0 and the conductivity as- 
sociated with the total electric field differs from that associated with the external 
field. Below we will calculate the conductivity and the conductance coefficients 
associated with the external field which may initially appear to yield a physically 
suspect result. However it must be remembered that the current induced by the 
external field already is linear in the field and can have no correction due to in- 
duced charge in linear response (see ref. [44]). So linear current response need 
only take into account correctly electron-electron interactions in the equilibrium 
state. For example the scattering cross-sections assigned to impurities must take 
into account self-consistent screening. However the currents are determined by 
the external potential and the voltages on the leads are fixed by the power sup- 
ply. Therefore despite the fact that standard texts often state that the conductivity 
associated with the total field is the only physically-relevant quantity [46], there 
is no need to calculate this conductivity in order to obtain the linear conductance 
coefficients. 

3.2. Kubo non-local conductivity tensor 

General many-body linear response theory [49] gives the induced current due to 
this perturbation as 

t 

J(r,t) = ^ liin / dt'e 61 ' cosujt' I dr' ( [j(i\ t), p(r\ t')]) 0(r'). 

k 6^0 J J 
— oo 

(3.2) 

Here the square brackets denote the commutator and the angle brackets de- 
note the thermal equilibrium quantum expectation values obtained by taking the 
trace of the operator multiplied by the density matrix of the grand canonical en- 
semble: exp[-/?K ]/^G where Zq — Tr{pc} is the grand partition function 
and K = H — fiN (H is the unperturbed hamiltonian including electron- 
electron interactions). The operators are written in the interaction representation 
and hence are time-dependent. Equation (3.2) expresses the current in terms of 
whereas to extract the conductivity tensor we wish to relate the current to the 
electric field. We can do this using a trick due to Kubo which is described in 
detail in ref. [46]. The expectation value of the commutator in Eq. (3.2) can 
be rearranged using the cyclic property of the trace and the fact that one can 
regard the operator exp[— /?K ] as generating displacements in imaginary time: 
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([j,Pg]) — _ P(t + The difference of the density operators is 

then expressed as the integral of p over an auxiliary variable A giving 

t (3 

J(r) = -elim J dt'e 6t ' cos tot' J dX j dv' ( j(r, t - ihX) p(r', t')) </>(r'), 

-oo 

(3.3) 

where we have used the cyclic property of the trace to shift the argument of j by 
— ihX. We can now apply the continuity equation ep(r, £') = —V • j(r, t') in the 
interaction picture to Eq. (3.3) and integrate by parts. This yields 

t & 

J(r,u;)=lim J dt' e 61 *' 1 cos ut' J dx(^(r,t - ihX) 



— oo 



X 



- J dr'i(v',t')- V'0(r') + J ds'j ± (s',t')(t>(s') 



(3.4) 



where jj_ is the normal component of the current density, and the surface s' 
must be taken to infinity because the original integral, Eq. (3.3), is of infinite 
range for the open system. In ref. [47] a careful argument is given demon- 
strating that the surface term vanishes as long as the correct order of limits is 
taken: lim,5_o linir— oo- Note that the bulk term is proportional to the elec- 
tric field as we desired and only extends over the finite sample region A where 
the field is non-zero. Since the remaining bulk term has the form J(r.u) — 
J drcr(a;, r, r') ■ E(r'), we can now identify the non-local conductivity tensor 
a(u) as an expectation value of the relevant product of current operators. This 
term can be expressed in terms of complete sets of (many-body) eigenstates 
|a), with energies E a ,Ep and the integrals over t' and A can then be done. 
Taking the limit to — > (but cot finite) we find a purely dissipative DC response 
(~ cos cot) with 

a(r,r ) = —in Inn > -r- — (3.5) 

v ' ; ^ E - E a E - E a + ihS 

a. j3 

where P a = (pc)aa is the thermodynamic occupation probability of many-body 
state a. As discussed above, this operator relates the current density to the external 
and not the total electric field, but this is sufficient for the purpose of obtaining 
the conductance coefficients. 

The non-local conductivity tensor cr(r, r') is a sample-specific quantity which 
may be defined for each impurity configuration, sample and measuring geometry. 
In general one does not expect this function to depend only on the difference r—r' 
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and one expects a spatially-varying current density J(r). If one averages this 
function over impurity configurations in the bulk one obtains of course a function 
of r — r' which is often taken to be short-ranged and of the form a(r, r') = 
cr exp[|r - r'\/l]\ thus the conductivity is taken to be local on a scale larger 
than the elastic mean free path. However this is not correct; Kane, Serota and 
Lee showed [53] that there are non-local contributions to both the average and 
variance of the Kubo tensor on the scale of the inelastic mean free path (leading 
to the non-local effects discussed in section (1.2)). 

It can be shown [40] that the sample-specific conductivity tensor of Eq. (3.5) 
satisfies the generalized Onsager relation cr;j(r, r' , B) — aji(r' . r, —B). It fol- 
lows that the spatially-averaged = J drdr'aij(r,r' , B) will satisfy the fa- 
miliar Onsager relation (by interchange of the dummy integration variables), but 
this spatially-averaged conductivity cannot be directly measured in mesoscopic 
samples. Hence the relevant time-reversal symmetry is the reciprocity symmetry 
discussed in section (2.3). 

3.3. Conductance coefficients in linear response 

The currents J m in Eq. (2.2) are obtained by integrating the current density ob- 
tained from Eqs. (3.4,3-5) through the cross-section of lead m. Once this is done, 
we must integrate again by parts with respect to r' in Eq. (3.4) to express the I m 
in terms of the voltages on the leads V m . The bulk term in this integration then 
involves J dS m V • a which can be shown to vanish from current conservation 
[40] (S m is the cross-sectional surface of lead ?n); the remaining surface terms 
give a contribution from each lead proportional to the voltage on that lead: 



where y m -Xm denote the transverse and longitudinal coordinates in lead m. 

We note that Eq. (3.6) involves a surface integral of exactly the same form as 
the one in Eq. (3.4) which we claimed vanishes, except that here there is no limit 
taking the surface to infinity. Obviously Eq. (3.6) must be non-vanishing for our 
linear response calculation to be sensible, which raises the question of whether 
our treatment of the earlier surface term was consistent. In fact Sols has recently 
claimed [50] that this treatment is inconsistent and that the continuity equation 
must be altered '"at infinity" to eliminate this earlier surface term. However the 
surface term in Eq. (3.4) required that the limit of the surface going to infinity be 
taken before the adiabatic limit 6 — > 0, whereas in (3.6) the limit S — + is taken 
and then (if desired) the limit of the surface tending to infinity may be taken. 
Following [47] we show explicitly below that these limits do not commute and 
that there is no inconsistency. 




(3.6) 
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Equations (3. 5), (3. 6) for the g mn are valid in the presence of arbitrary interac- 
tions, but are not so useful microscopically in the general case. Thus we special- 
ize to the case of non-interacting quasi-particles. In the non-interacting limit, the 
statistical weights P a can be replaced by Fermi functions f{E a ), and |a) become 
single-particle states. Once this is done equation (3.5) is often divided into two 
terms by applying the Cauchy principal value identity to the energy denominator: 
[E/3 -E a + ikd}- 1 = -Trih6(E — E a ) + P{[Ep — Ea] -1 }. The delta function 
term is symmetric in magnetic field (without interchanging r — > r') while the 
principal value term is anti-symmetric, hence it must vanish at B = [40]. The 
delta function term only involves states within kT of the fermi surface, whereas 
the principal value term appears to depend on all states beneath the fermi surface 
(and it is sometimes stated that the Hall conductivity arising from this term is 
not a fermi-surface quantity). By deriving Eq. (2.2) below we will demonstrate 
that the Hall resistance is in fact completely expressible in terms of transmission 
coefficients at the fermi surface, a non-trivial result. To achieve this it turns out 
to be more convenient to work with the sum of the symmetric and anti-symmetric 
parts as in Eq. (3.5). 

The relevant single-particle states for our model are scattering states which may 
be expressed in the leads in terms of the elements of the scattering matrix S; in 
lead /, they have the asymptotic form 

= 1>Eap = *Vl ta + E S W« ■ W 

a' 

Here, the label a consists of the energy E, subband index a and lead p of the 
incident wave. An analogous definition holds for ip l p = ift l E , bq . The quantum 
wire eigenfunctions have the form 



4a(r) 



dE. 



-1/2 



e tk W± a (y), (3-8) 



where x± a are me transverse wavefunctions and k z ±a is the outgoing (+) or in- 
coming (-) wavenumber. With the above choice of normalization, the symbolic 
sums over collective labels have the explicit form Y^ a ~ * S P D a / ^ n me 
following, subscripts a imply a dependence on the energy variable E, while in- 
dices b belong to energies E' . As was shown in Ref. [40], current conservation 
implies the following properties for two wire eigenfunctions at the same energy 
in the asymptotic region: 

/ dm (& a \jM\Ca') = ^<W<W (3-9) 

with the notation u — ±, and round brackets for the single particle current matrix 
elements. If e l n is the propagation threshold of subband a in lead /, the summations 
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Jk lm(E' ) 



«■ 



*- 



Re(E ' ) 
-» > 



Fig. 6. The two possible contours that yield exponential decay in x' off the real axis for the various 
terms in Eq. (3.13). Open dots: two possible locations of the pole at E' — E — ih.6. 



in Eq. (3.5) lead to integrals of the type 

,f{E')-f{E) F S (E\E) 




E' - E E' — E + ihS ' 
where 



(3.10) 



F S (E\E) = J d*' (^ fcg | j j. (r) |^ ap )(^a„ b'x (s') I^S'fc, ) (3-11) 

contains the current matrix elements and thus vanishes at threshold. The surface 
S of integration is the surface of the sample region, 9 A Although we wish to let 
dA — > oo eventually, we can do so only after 8 — ► 0. 

We now evaluate the integral over E' in Eq. (3.10) by Cauchy's theorem; a 
similar approach is taken in Refs. [51,52]. In general the integrand in Eq. (3.10) 
may have a complicated singularity structure in the complex plane, arising both 
from the "Matsubara" poles of the first factor and the singularities of the S-matrix 
contained in Fs(E,E'). However since asymptotically the wavefunctions have 
a plane-wave dependence on x' (the longitudinal coordinate in the lead), we shall 
see that contributions from the return contour may always be made to vanish as 
long as the correct half-plane is chosen to close the contour. 

Since the locations of all the singularities of the integrand in Eq. (3.10) are 
independent of both the asymptotic and adiabatic limits except for the pole at 
E' = E — ihS, for convenience we choose contours C\, (shown schematically 
in Fig. ) which by assumption enclose only the latter pole. Hence C\ will contain 
no singularities at all, and terms which in the asymptotic limit may be closed 
in the upper-half plane will give zero contribution. Co will enclose the pole at 
E' — E — ihS unless it approaches the real axis outside the range of the E' 
integration (i.e. unless E is less than the sub-band threshold relevant for the term 
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in question). By this convenient choice of contours we need only evaluate the 
residue at E' — E — ihb. 

The residue theorem yields with v — 1 , 2 



dE + J dE ) E'-E E'-E + il 



'I 



6. ' - - - ~ (3.12) 

-2xi«,. a 0(e - 4) M-w-m Fs{E _ m E)t 

—i no 

where the step function 6 enters since C2 encloses a pole only if E > e b , and 8 U 2 
reflects the fact that no poles are enclosed by C\. Along C v , the limit S —> 00 
enables us to apply the asymptotic forms of ip a and ipp in the second current 
matrix element in Fs(E',E) to get in lead I 

(ta\j±\t b )+S p l^S lq Mta\jA.\e +b >) 

b' 

a' a'b' 

(3.13) 

irrespective of whether 6 = or not in Eq. (3.12). For E' in the complex plane, 
the imaginary parts of the wavenumbers satisfy 

sgn (lmk l ±b (E')) = ±sgn(ImE') . (3.14) 

This can easily be checked for the special case k l ±b oc ±y/E' -e q b , but holds even 
in asymmetric leads at B ^ where E' = e q h can occur for nonzero k. Note that 
the branch point at e q h causes no problems since Fs(E' ,E) vanishes there. The 

first and third term of Eq. (3.13) depend on £ l _ b oc e lk -' ,x , so that they acquire 
an exponential decay factor e - K ( E ') x ' when integrated along the contour C2. For 
the other two terms in Eq. (3.13), we choose C\ to obtain positive imaginary parts 
in k l +b (E'). The asymptotic limit x' — » 00 causes the contributions from C\ and 
Co to vanish identically. The result is that the real axis integrals over E' in Eq. 
(3.10) vanish due to Cauchy's theorem unless C2 encloses a pole. Its residue 
in Eq. (3.12) contains a factor e^Ut^-^)*'. For E / e q b we can expand the 
outgoing wavenumber as k l _ b (E) — i8 v~ l (E) where the group velocity v(E) — 
(l/h)dE/dk is negative. We thus pick up a decaying exponential of the form 

F S (E - ih6,E) oc exp [-6 x'] . (3.15) 

This exhibits clearly the noncommutativity of the limiting procedures: If we take 
x' — ► 00 first as we were required to do in the surface term of Eq. (3.4) then this 
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term vanishes. On the other hand, the g mn of Eq. (3.6) are nonzero because there 
we are required to take <5 — > first in Eq. (3.15), yielding 

oo 

lim , dE , /(£')" HE) F S (E>,E) 



s~o J E'-E E'-E + ih6 (3.16) 

= "-2m 6 Vt2 Q(E - el) f'(E) E S (E, E). 

This result leads to Eq. (2.2), as we now show by taking the asymptotic limit 
5* — ► co. The product of current matrix elements at the same energy in Fs{E. E) 
can be evaluated with Eq. (3.13) if we let r go to infinity as well. Considering 
only n ^ m. 7 of the 16 terms in F${E, E) vanish identically due to Kronecker 
delta factors like 8 pm 6 pn . When the sums over the incident wave parameters p, a 
or q, b are performed five more terms drop out due to the unitarity of the S matrix 
[47], two of the remaining four vanish due to integration along Ci (see 3.16) and 
one vanishes due to Eq. (3.9). 

Thus our expression for g mn now contains only one term: 



N oo 



oo 



gmn = ih j dy m J dy n ]T ]T (-2™) J dEQ(E - e q b ) f'{E) 

p, ( } = l a,b=\ v 
X^gn<5 p n^ [S*nn,b' bSmn,a' a (£+6' \j± \€+ a ' ) (£-a \j ± \€-b) ] £> = l 



e 2 



a'b 1 

oo 



= - / dE [-/'(£)]£ \S mn , a - a \\ (3.17) 



where Eq. (3.9) has been used. This is the desired result, Eq. (2.2). We note that 
the term identical to this one with the lead indices reversed vanished above due 
to the sign of 6 (which required integration along C\)\ i.e. due to the adiabatic 
switching on of the perturbation from t = — oo. So the causality requirement that 
g mn depend only on scattering probabilities from n to m, and not vice versa is 
correctly represented by this procedure 



4. Ballistic transport phenomena 

4.1. Novel phenomena in high-mobility microstructures 

We have shown above using physical arguments and linear response theory that 
the transport properties of mesoscopic conductors can be plausibly related to a 
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quantum many-channel scattering problem. One of the first and most import- 
ant applications of this scattering approach was to the Aharonov-Bohm effect 
[54,55] and the universal conductance fluctuations in disordered microstructures 
[8,9,13,19]. These mesoscopic interference phenomena in disordered systems 
are reviewed in the chapters by Altshuler [16] and Bouchiat [15]. For disordered 
systems the analytic technique of impurity-averaged perturbation theory is ap- 
plicable and leads to powerful results (such as the universality of the conduct- 
ance fluctuations). The scattering approach is primarily useful here for numerical 
simulations. On the other hand, starting roughly in 1987 [31,32,56,41,57,43,58- 
61], reduced dimension high-mobility GaAs microstructures became available for 
experimental study and a wide variety of novel transport phenomena were dis- 
covered. Since measurements could be made in such structures on a scale shorter 
than the bulk elastic mean free path, electronic motion was ballistic (at least over 
several collisions with the boundaries). In this case the scattering approach was 
found to be particularly useful since it allowed a simple modelling and in many 
cases a simple physical understanding of of the novel phenomena. We will review 
several examples in detail below. 

In such microstructures the potential confining the electrons is approximately 
uniform in each straight segment of the structure so that the concepts of sub-bands 
and propagating channels apply within the structure itself (and not just in some 
idealized perfect leads). For this reason such segments are sometimes referred 
to as "quantum wires" or "electron waveguides". We have already discussed in 
section (2.4) the simplest consequence of such sub-band structure, the quantized 
two-probe point contact resistance [31,32]. Although multi-probe structures do 
not show resistance steps associated with sub-band threshholds, they do show 
sharp resistance minima at low temperature [61]. These fermi energy dependent 
features represent threshhold effects in the quantum S-matrix of the type already 
well-known in nuclear physics. Due to lack of space they will not be treated in 
any detail below. 

The type of samples we shall consider in this section are GaAs microstruc- 
tures with quasi-one-dimensional leads which either intersect to form a four- 
probe junction [56,58] (see Fig. (7)) or which enter and leave a larger 2D is- 
land (see insets to Figs. (10),(12)) of electrons [62,64,65]. In the case in which 
the island is strongly coupled to the leads it acts as a resonant cavity for the 
transmitted and reflected electrons; when it is isolated from the leads by a tun- 
nel barrier it acts in some respects like an artificial atom, often referred to as a 
"quantum dot". We shall refer to the four-probe junctions and two-probe res- 
onators as quantum junctions, reserving the term quantum dot for the tunnel- 
ing systems. The quantum junctions are simpler in that straightforward use of 
the scattering approach for non-interacting electrons appears to describe many of 
their properties adequately. In the quantum dots the charging energy [62,63,66] 
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and Coulomb blockade [67] phenomenon determine the basic experimental ob- 
servations. We shall not have space to treat the quantum dots in any detail, al- 
though we have argued that quantum interference of single electrons still plays a 
role [68]. 

The experimental transport phenomena discovered in the ballistic quantum 
junctions can be roughly divided into three categories: 1) Classical Geometric 
Effects which can be understood from classical scattering of ballistic electrons 
from the boundaries or geometric features of the junction. 2) Quantum Thresh- 
hold Effects which we will not discuss further. 3) Quantum Interference Effects 
which arise from interference of different electronic scattering trajectories. This 
interference gives rise to both conductance fluctuations (CF) and weak localiza- 
tion (WL) even though the conductor has negligible bulk disorder. 

Although in the most general sense the ballistic WL and CF effects have the 
same physical origin as in disordered systems (which may be regarded as chaotic 
systems specified by a large number of parameters), these interference effects 
are interesting from the perspective of the theory of quantum chaos (quantum 
manifestations of chaotic classical dynamics) [69]. There are very few quantum 
systems which can be studied experimentally and which have simple enough dy- 
namics that a detailed connection can be made between the classical and quantum 
properties via semiclassical quantum theory. Jalabert et al. [70] were the first to 
point out that ballistic quantum junctions might provide such a class of systems. 
In this case the statistical properties of the conductance were predicted to reflect 
the nature of the classical mechanics. If the system is fully chaotic certain statist- 
ical behavior has been predicted both for the WL and CF which will be discussed 
in detail below. But perhaps even more intriguing was the possibility of determ- 
ining the nature of the classical dynamics from the analysis of the conductance 
statistics; i.e. distinguishing chaotic from regular dynamics. Since the underlying 
electron system has substantial electron-electron interactions, such experimental 
results would provide a strong confirmation of the quasi-particle picture of fermi 
liquid theory, in a regime in which it has never before been tested. In fact we 
shall see below that there is now fairly solid experimental support for the idea that 
mesoscopic conductance can distinguish classically chaotic from regular electron 
dynamics. In the next section we will develop the semiclassical formalism which 
can be used to understand both the classical and quantum-chaotic effects in these 
systems. 

4.2. Semiclassical formalism 

Our aim is to combine semiclassical quantum theory with the scattering approach 
of Eq. (2.2) to develop a semiclassical theory of ballistic junctions. Thus the 
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basic objects of interest are the conductance coefficients: 





(4.1) 



/ y -^mn,afe) 
a, 6=1 



where a, b label the outgoing and incident channels in the leads. We will sup- 
press the lead indices m, n henceforth and focus on the two-probe case (where 
the transmission coefficient directly yields the conductance) for simplicity. As- 
sume a two-probe quantum juntion with the left lead joining the junction along 
the line x = and the right lead joining along x = L(left inset, Fig. (12)). The 
transmission amplitude from a mode b on the left to a mode a on the right for 
electrons at the Fermi energy E is given by the projection of the retarded Green 
function (evaluated at the Fermi energy) over the transverse wavefunctions 4>b 
and <f) a of the incoming and outgoing modes [28]: 



t ab = -ih{v a v h ) 1 ' 2 / dy' / dy 4>* a {y') (f> b {y) G(L, y'; 0, y; E). (4.2) 



Here the integrals over y, y' are over the transverse cross-sections of the left and 
right leads as in section (3.3); the reflection amplitude is of the same form except 
that integrals are over the same lead and there is an additional term <5 a & represent- 
ing the incoming wave. 

This equation is exact and is the basis for the numerical results given below 
in which the exact Green function is found using a recursive algorithm which 
is described in detail in ref. [71]. It also provides an appropriate starting point 
for the semiclassical approximation to the transmission and reflection amplitudes. 
We proceed by replacing the Green function G in Eq. (4.2) by its semiclassical 
path- integral expression [69] 



given as a sum over classical trajectories s between points y and y' of the entrance 
and exit cross sections (hereafter we will not write any x dependence). Here 6 
and 0' are the incoming and outgoing angles, D s = (m/v\ cos#'|) | (d9/dy') y |, 
S s is the action integral along the path, and /i is the Maslov [69] index given by 
the number of constant-energy conjugate points. 

In the case of hard-wall leads, the transverse wavefunctions have the form: 
<t>b{y) = \/2/W sm(biry/W). For large integers b the integral over y in Eq. 
(4.2) will be dominated by the stationary-phase contribution occurring for traject- 




G sc] (y' ]y] E) = 



(2tt^)( 3 / 2 ) 
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ories starting at points yo defined by 
fdS\ bh-n T 

{dy-) y r- p » = -w> 6 = ±6 ' (4 - 4) 

that is, the dominant trajectories are those for which the initial transverse mo- 
mentum equals the momentum of the transverse wavefunction. Performing also 
the y' integral by stationary-phase, we obtain 



tab = s S n (^) s S n (^) ft** exp (-S s {a,b,E) - i-ft^\ , 

s(a,b) ^ - / 

(4.5) 

where the sum is now taken over trajectories s between the entrance and exit cross 
sections with incoming and outgoing angles 9 and 9' such that sin 9 = Jm/kW 
and sin 9' — air/kW. The reduced action is 

5(5, 5, E) = S(y' ,y . E F ) + hnby /W - fmay^/W, (4.6) 

the new pre-exponential factor is 



b.= 1 



mv\ cos 9' 



dy_ 
d9' 



(4.7) 



(The explicit form of the new Maslov index, fi, is known [72] but will not be 
needed here.) In general there will exist families of trajectories which directly 
traverse the junction without scattering from the boundary for which the second 
stationary phase integration is invalid [72]; these direct trajectories lead to non- 
universal effects which we neglect in our discussion below. Experiments and 
simulations have typically employed geometries intended to minimize these ef- 
fects. 

Equations (4.1) and (4.5) express the conductance as a sum over certain pairs 
of classical paths. If we now specialize to junctions with a constant potential in- 
side and hard wall boundaries ("billiards") the classical trajectories do not change 
with energy (or wavevector k — \f2mE/Ti) and it is convenient to express the 
semiclassical conductance as 

T (k) = £^ = ~EEE cr<*>. < 4 - 8 > 

a.b= [ a,b s u 



K^ k ) = >/IXexp [ik(Ls ~ In) + Z7T0 S , U ], (4.9) 

where s and u label the paths with entrance and exit angles sin# = bn/kW 
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and sinfl' = air/kW, the effective length is L s = S s /kh = S s /kh + ysind- 
y' sin 9', the Maslov indices are included in the phase-factor SiU (which we sup- 
press henceforth), and A s = (hk/W)D s . Note that by these definitions A is 
independent of energy so that the only energy-dependence in the summands is 
the explicit k-dependence in the actions. 



4.3. Classical two-probe conductance 



The classical conductance may be obtained from the semiclassical formula (4.8) 
by approximating the mode summations by integrals over continuous angles 6, 6'. 
That this is correct may be seen by noting that the momentum p = hk is a classical 
quantity which must remain fixed in the classical limit (/?, — > 0), thus requiring 
the limit k — > oo at the same time. Since the entrance angle sin# = bn/kW, 
the allowed angles become continuous in the classical limit. Noting further 
that the number of channel N c = Int{fcW/7r} we arrive at the prescription 
(rc/kW) ~~ *■ S-\ dsin# and similarly for sm9' (the lower limit is —1 be- 
cause the summation in Eq. (4.5) is over positive and negative a, b). With this 
substitution Eq. (4.8) becomes 




dsm9 / dO' V 

-1 J-n/2 




1/2 

% ~ ~ 

exp [-(5 S - S u )}. 

(4.10) 



Equation (4.1 1) indicates that some further thought is required to get a mean- 
ingful classical limit for the conductance. Obviously linifc-^ T(k) is not well- 
defined since it diverges as k and apparently oscillates indefinitely due to the 
k-dependent phases in the summation. The divergence is real in the sense that 
the conductance of the ballistic microstructure really is proportional to N c = 
Int{kW /n}. The useful classical quantity to describe the system is not the trans- 
mission coefficient (which is normalized to N c ), but the transmission probability 
normalized to unity for perfect transmission T = T/(kW/ir). This quantity is 
finite as k — > oo but still will have the rapid phase oscillations for any finite k. 
These rapid oscillations (which are just the conductance fluctuations with energy) 
will turn out to be order unity in T implying that they are order l/N c in T; hence 
they will tend to zero as k — > oo. However a convenient and physical method 
of eliminating them at this point is to define the average of T over wavevector 
k. This method is physical because a non-zero temperature will have essentially 
the same effect and the classical ballistic effects of interest are observed pre- 
cisely at intermediate temperatures for which the inelastic scattering length is of 
order the size of the microstructure but the temperature is high enough to aver- 
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age out the quantum interference effects. The average of a quantity A(k) will be 
defined as: 

(A)= lim(l/<7) r C + ' dkA(k), q c W/7c^>\. (4.11) 

q ~°° Jq c 

Note that in principle we begin averaging at q c > n/W so that the semiclassical 
approximation is valid over the entire range. If we divide Eq. (4.10) by kW/n 
and perform the k average then since k only appears in the exponents we obtain 
rigorously [72] zero contribution unless the paths s, u have the same actions. If 
we neglect for the moment paths with the same action due to an exact symmetry 
(to be discussed in the WL section below) then only terms with s = u contribute 
to the sum. Noting that for the diagonal terms s = u the integral over outgoing 
angle 9' can be combined with the factor dy s /d9' to yield an integral over the 
initial position y, the remaining integrals can be expressed in the form: 

T = \ d(sin 60 jf (4 - 12) 

where /(y. 9) = 1 if the trajectory with initial conditions (y, 9) is transmitted and 
f(y,9) = otherwise. T in Eq. (4.12) may be calculated by a completely clas- 
sical simulation, i.e. by injecting an ensemble of point masses into the junction 
and counting the fraction transmitted for a given magnetic field and geometry. 
The "classical" approximation to the quantum conductance is just 

e 2 kWT 

g=-, • (4.13) 

ll 7T 

This equation generalizes to the multi-probe conductance matrix g mn simply by 
the replacement T — > T mn ; it was first introduced to describe ballistic junctions 
by Beenakker and van Houten [73] using a different argument. Equation (4.13) 
may also be derived from a Boltzmann equation approach [71]. The notion that 
the classical transmission probability can be used in the two-probe or multi-probe 
conductance formulas is not particularly surprising. However the present argu- 
ment allows one to determine the correct initial distribution of injected particles 
(uniform in space and with weight cos<9 in momentum) and more importantly 
provides a natural generalization to the WL and CF effects. 

4.4. Classical effects in ballistic junctions 

Perhaps the first novel transport phenomenon discovered experimentally in bal- 
listic junctions was "quenching" of the low-field Hall resistance. Roukes et al. 
[56] and subsequently several other groups [57-59] observed that the low field 
Hall resistance of ballistic junctions was strongly suppressed up to some field at 
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which it rose sharply, overshooting the bulk 2D value and flattening out again 
to form a plateau (similar to quantum Hall plateaus but not at a quantized value 
of Rh). This effect as measured in the later experiment by Ford et al. [58] is 
shown in Fig. 4. Although initially a number of theoretical explanations for the 
effect were proposed, it is now fairly well accepted that it is due to geometric 
features of the junction and has an essentially classical origin (in the sense that 
Eq. (4.13) is classical). The particular property of the junction which was found 
to lead to quenching was the widening of the leads as the junction is approached 
(see Fig. 7). This widening is almost unavoidable experimentally because the 
conducting regions of the GaAs microstructure are defined by electrostatic gates 
placed in layers substantially above the two-dimensional electron gas. The spacer 
layer between the 2DEG and the gates is essential to get the high mobility, but 
it means that geometric features patterned in the gate are smoothed out at the 
level of the 2DEG. The widening of the leads creates a tendency for the entering 
electrons to have their momentum rotated from the transverse to the longitudinal 
direction (see Fig. 7), i.e. the electrons are collimated in the forward direction 
as they enter the junction. For this to occur the widening must be "adiabatic": 
the leads must widen gradually enough that the change in transverse momentum 
between collisions with the side wall is small compared to the initial transverse 
momentum [36,75,76]. Baranger and Stone [74] were the first to point out that 
this collimation could explain the quenching phenomenon. 

If one assumes that the four-probe ballistic junction has four-fold rotational 
symmetry (as was approximately true in the experimental structures) the Hall 
resistance of Eq. (2.2) can be written in a simple form in terms of the transmission 
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coefficients Tr,Tl,Tf (transmission coefficients for the right,left and forward 
lead): 

P h (T R + T L )(T R -T L ) 

R H = ~j n(rr rp rp V • (4.14) 

Here the denominator D{Tr,Ti,Tf) is a sub-determinant of the conductance 
matrix [74]; it is relatively insensitive to magnetic field and its explicit form will 
not be needed. Note that Rh is proportional to the asymmetry Tr—Tl in left-right 
transmission induced by the field, as one might expect. However Rh is also pro- 
portional to the total sideways transmission Tr 4- Ti\ Baranger and Stone argued 
that this factor is suppressed by the collimation effect induced by the widening of 
the leads. They were able to establish [74] by quantum simulations of equation 
(4. 14) that widened leads were essential to obtain generic quenching (for example 
a perfect square cross junction did not show quenching after averaging over en- 
ergies (see Fig. 8)). The simulations at T — also showed the large quantum 
conductance fluctuations which will be the subject of section (4.9) and it was only 
after including the effects of temperature by energy-averaging that generically flat 
low field Hall resistance was found (again see Fig. 8). 

The simulations also clearly indicated that quenching did not arise solely from 
the suppression of sideways transmission [74]; it was found numerically that the 
widening of the leads was also reducing the asymmetry factor in Rh, Tr — Ti. 
At roughly the same time and independently, Ford et al. [58] proposed a simple 
mechanism for the asymmetry reduction. If the junction widens gradually elec- 
trons diverted by the field towards the left probe might rebound into the opposite 
("wrong") probe (see inset to Fig. 4). The experimental data of Fig. 4, which 
were obtained from junctions intentionally fabricated to enhance this effect, dra- 
matically confirm it. At an intermediate field the sign of the asymmetry is actually 
reversed, leading to an inverted Hall resistance (i.e. voltage induced opposite to 
the Lorentz force). This and other related experiments [58,59] established that 
junction geometry was crucial. 

Shortly thereafter Beenakker and van Houten [73] proposed that the quench- 
ing, inversion and last plateau observed in the Hall resistance of ballistic junctions 
could all be explained in terms of Eq. (4. 13) as classical magneto-size effects. By 
classical simulations of the type described after Eq. (4.12) they were able to re- 
produce the experimental results qualitatively and semiquantitatively. They went 
on to suggest that yet a third mechanism might be important in the quenching 
which they referred to as "scrambling". The idea was that the rounded junction al- 
lowed long trapped trajectories which would "forget" which lead they had entered 
from and hence which Hall probe the field was favoring. Subsequently Baranger 
et al. [71] performed detailed classical and quantum simulations of four-probe 
junctions looking at both the Hall and longitudinal resistance (which is often re- 
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ferred to as the "bend" resistance since the current turns a corner). They looked 
carefully at the three classical mechanisms for quenching: 1) Forward transmis- 
sion enhancement. 2) Rebound mechanism. 3) Scrambling. It was shown from 
classical simulations that the asymmetry Tr — Ti was determined almost com- 
pletely by short trajectories and that while the rebound and forward enhancement 
effects were both important for the quenching, the scrambling mechanism was 
not significant [77]. The classical model for the last plateau and bend resistance 
also achieved reasonable success. In summary, the classical conductance formu- 
las Eqs. (4.12), (4.13) provide a simple and reasonably good description of the 
magnetotransport through ballistic junctions, in the temperature range 1 — 4AT. 

4.5. Quantum ballistic effects 

The classical theory works well at intermediate temperatures but begins to fail 
badly as the temperature is lowered another order of magnitude. One source of 
the failure is the quantum threshhold effects [61] mentioned above which become 
more pronounced as temperature smearing is reduced. However the more inter- 
esting source of failure is the appearance of interference effects associated with 
multiple bounce trajectories trapped in the junction. In Fig. 8 (top) we show 
the Hall resistance of a ballistic junction at lOOmK; although one can make out a 
tendency towards quenching, superimposed on this smooth classical magneto-size 
effect are reproducible fluctuations similar to the UCF of disordered microstruc- 
tures. The presence of the quenched envelope in the data is important because 
it indicates that the transport is truly ballistic and hence the fluctuations are not 
simply due to disorder arising from poor sample fabrication. In Fig. 8 (bottom) 
we show theoretical calculations which compare the classical and quantum pre- 
dictions for the the Hall resistance of ballistic junctions. As suggested by the 
experiment of Fig. 8, the classical simulation reproduces the quenching behavior 
but not the fluctuations, which in the quantum simulations occur both as a function 
of magnetic field and energy. Simulations demonstrate that a rounded four-probe 
junction definitely generates chaotic scattering for the trajectories which remain in 
the junction for many bounces [70,78]. But the classical transport coefficients of 
Eq. (4.12) involve a phase-space average over many trajectories which smooths 
out the strong dependence on external parameters of the individual chaotic traject- 
ories. Classical chaos is not easily observed in ballistic microstructures. Although 
the fluctuations were observed experimentally as early as 1988 [57], it wasn't un- 
til after their connection with quantum chaos was proposed in 1990 [70] that they 
became the object of intensive study [64,65,79-83]. However the occurence of 
the fluctuations is not a unique signature of quantum chaos. Both chaotic and in- 
tegrate cavities exhibit such fluctuations [72]; it is only their statistical properties 
which differ in the two cases (as we shall see in detail below). 
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B/B 

Fig. 8. (Top) Experimental Hall resistance from ballistic junction at 100 mK (adapted from [60]). 
(Bottom) Four calculations of the T = Hall resistance of a ballistic junction: (dashed) Classical 
calculation from Eq. (4.12) showing quench but no fluctuations, (dash-dot) Quantum square-corner 
junction (R/W = 0) from Eq. (4.1) showing no fluctuations due to lack of trapped trajectories. 
(Solid line) Quantum rounded junction shov/ing fluctuations. (Squares) Energy-averaged quantum 
calculations for the rounded junctions showing quenching. Inset: Schematic of rounded junction, 
results are for R/W = 4. /?.// is in units /i/e 2 , B in units Bq = mcvj/eW. 



Keller et al. [65] were the first to fabricate microstructures of the relevant type 
in which the electron density (and hence the fermi momentum kj) was tunable 
while maintaining the geometry approximately fixed. Some representative data 
from their measurements is shown in Fig. 9. Here ballistic conductance fluctu- 
ations are plotted vs. k / for several values of the magnetic field. One sees that 
on average increasing the B field increases the conductance: this is the ballistic 
weak localization effect. It is more difficult to observe in ballistic microstruc- 
tures because the averaging out of the fluctuations must be done explicitly (either 
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Fig. 9. Measured conductance fluctuations as a function of fey in a ballistic junction at magnetic fields 
varying from 35 to 350 G. AG denotes the difference G(B, kf) — G(0, fc/) and is typically positive 
due to WL effect (Adapted from [83]). 

by varying the fermi energy [65] or by making many similar samples in paral- 
lel [81]), in contrast to disordered systems where studying samples much longer 
than the phase coherence length leads to self-averaging of the fluctuations. Figs. 
8, 9 shows that both the WL and CF effects exist in the ballistic regime. As in 
disordered microstructures the origin of the effects is constructive interference of 
time-reversed pairs of electron trajectories in the case of WL [4,16,84] and in- 
terference of more general trajectories with a field and energy-dependent phase 
for the case of CF [8]. However in the ballistic case a more detailed semiclas- 
sical theory is possible. Also it is now understood that there are subtle statistical 
differences between the disordered and ballistic effects (see section 4.11 below). 

4.6. Ballistic weak localization: coherent backscattering 

Although the conductance fluctuations are the most prominent feature of the data 
in Figs. (8),(9), the weak localization effect is a simpler starting point for the 
semiclassical theory and demonstrates both its strengths and weaknesses. In the 
discussion of the classical conductance in section (4.3) the interference of unlike 
paths in Eq. (4.10) was eliminated by the wavevector average (Eq. (4.1 1)) which 
introduced a delta function setting the effective actions S s — S u . At that point we 
neglected the possibility that unlike paths may have the same effective actions. 
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This may occur due to exact spatial symmetries of the confining geometry (we 
have discussed this case elsewhere [71]) but more importantly it will occur due to 
time-reversal symmetry. The breaking of this symmetry by a weak magnetic field 
leads to the WL effect in conductance which we now attempt to evaluate using 
the semiclassical theory. 

In doing so we immediately encounter a conceptual difficulty. The semiclas- 
sical transmission amplitude t a t, of Eq. (4.5) is given by a sum over classical paths 
entering at the quantized angle sin 6b = ±bn/W and exiting at sin 6' a = ±air/W; 
the time-reverse of such a path exits at sin 6' b and enters at sin 6 a . Hence the time- 
reversed path contributes to tf, a and does not interfere with the original path in 
the summation of Eq. (4.8) (which is over |£ a &| 2 ). The only S-matrix elements 
in which exact time-reversed pairs interfere semiclassically are the diagonal ele- 
ments of the reflection matrix, R aa = |7' aa | 2 . Hence it is only the effect of TR 
symmetry on Rd = ^2^ ' mat we can evaluate using the method of section 
(4.3). Obviously Rd is related to the transmission and off-diagonal reflection 
(Rod — J2a±b R"t>) by unitarity: Rd + Rod + T = N c . Thus the failure of our 
semiclassical approach to find a direct WL effect in transmission is troubling; we 
shall return to this point below. 

We first note that the semiclassical reflection amplitudes have exactly the same 
form as Eq. (4.5) except that both and & refer to the same lead and there is 
a term 8 a b representing the incoming wave. Since below we will be concerned 
only with the field-dependent correction to Rd we can neglect the 6 a t term. The 
relevant part of Ro{k) in semiclassical approximation is: 

SD = i>„a = i^££F*,vm (4.15) 

a a s,u 

where F*-%(k) is given by Eq. (4.9) generalized for reflection. Rd has a con- 
tribution from like paths s = u of the type leading to Eq. (4.12); however this 
contribution is not sensitive to weak fields. Hence we immediately focus on the 
terms s — T(u) in Eq. (4.15)(where T is the time-reversal operator). These 
terms will survive the wavevector average and give an order unity contribution 
to (R(k)) which we denote by SR D {B). At B=0 S T ( U ) = Su, Ar(u) = D u and 
since | sin 6\ = | sin 6'\ for paths contributing to R aa we can again use the relation 
(n/kW) ^ c d sin 6 to obtain 

6R D (B = 0) = ^ j d{smd) ]T A s = i / d(sin 6)A{0, 6). 

1 s(0, ±6) ~ J ~ l 

(4.16) 

The quantity A(6,0) - J2 s (e ±0)(|W COS measures the stability 
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of classical paths which enter and exit at 9 with respect to variations in the ini- 
tial position. Note that in the calculation for T (the terms s = u in Eq. (4.10)) 
the same quantity appears generalized to unequal #, 9' , A(9, 9'). For chaotic clas- 
sical scattering there is no reason to expect A(9, 9') to depend strongly on its argu- 
ments. Hence we make a uniformity hypothesis: A(9, 9) « d(s'm 9')A(9, 9'). 
Classical simulations indicate that this assumption is reasonably well satisfied in 
chaotic cavities without direct trajectories [72]. If we substitute the uniformity 
relation into Eq. (4.16) it becomes identical to 7Z, the classical reflection prob- 
ability (defined in obvious analogy to T of Eq. (4.12)). Since the terms s = u 
in the semiclassical sum will give exactly the same contribution to Rp at B = 
this implies 

R D (B = 0) « 211. (4.17) 

Since off-diagonal reflection coefficients and transmission coefficients will not 
have the contribution from time-reversed pairs semiclassically we conclude that 
the diagonal reflection coefficients are on average twice as large as these other 
elements in the S-matrix. This factor of two enhancement of the backscattering 
due to phase coherence and time-reversal symmetry (also known as the elastic 
enhancement factor in nuclear physics) is well-known in the literature on weak 
localization [4]. 

4.7. Ballistic weak localization: field dependence 

We now evaluate the effect of a magnetic field on the enhanced diagonal reflec- 
tion. Introducing a magnetic field in general changes both the classical paths 
traversed and the action along a given path; however, we will be considering 
low-field effects in which the change in the geometry of the paths is negligible 
and only the phase difference which now appears between time-reversed paths is 
important. For time-reversed pairs this phase difference (which is essentially due 
to the Aharonov-Bohm effect) arises due to the different sign of the "enclosed 
flux", (S 3 - S u )/h = 20 S £/0 O where S = (2tt/£) f s A-dl is the effective 
"area" enclosed by the path and O = hc/e. Although we refer to S as an area, 
the paths we are considering are not in general closed loops and thus S is not 
itself a gauge-invariant quantity. We have checked that the statistical properties 
of Q s are gauge-invariant in the chaotic case. Neglecting the change in classical 
paths we find instead of Eq. (4.16) 

6R D (B) = \ f d(sin0) £ A s e l2e ^ B ^\ (4.18) 

s(0, ±8) 

We now wish to convert the sum over classical paths s to an integral over the 
effective area of the paths 0. Note that in the summation a classical quantity A s 
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is multiplied by a phase factor with a quantum scale O = hc/e so as h — > 
A s may be treated as slowly-varying with respect to the phase factor when we 
convert from a sum to an integral. Hence we obtain 

/OO 
d0P(0)exp[i20B/0 o ], (4.19) 
-oo 

where P(0) is the probability that a traversing classical path encloses an effective 
area 0. 

Both numerical calculations [70] and analytic arguments [85,86] have found 
that the area distribution for long orbits takes on a universal form for simple 
chaotic billiards: P(0) oc exp (-o. c /|0|), where the only parameter character- 
izing the classical phase space is the quantity a c ;, the inverse of the typical area 
enclosed by a scattering trajectory. The origin of this universal form may be un- 
derstood qualitatively by the following argument [85,86,72]. First, the number 
of orbits which remain in the scattering region long enough to cover a distance 
L decays exponentially with L. Crudely the reason is that for long chaotic or- 
bits the probability of escaping at each encounter with the boundary is a constant 
p <C 1 given approximately by the ratio of the size of the opening to the peri- 
meter of the billiard and is uncorrected with previous bounces. Hence after n 
bounces P(n) ~ e -nlnp , and if is the mean length traversed between bounces, 
P(I) ~ e -(L/d)\n P = e -idL Secondi the 

area Q s of a given trajectory is es- 
sentially given by its winding number around a central point in the billiard times 
a typical area per circulation A (it is precisely this if the field is replaced by an 
Aharonov-Bohm flux, as considered by Berry and Robnik [87]). For a chaotic 
system the area accumulates diffusively, as both senses of circulation occur with 
roughly equal probability for long orbits. Hence the mean winding number for an 
orbit of n bounces is zero and its variance w\ ~ n [87,85,86]. The distribution 
of areas for orbits of length L is then Gaussian with the variance A 2 « {L/d)A^. 
Finally, in order to obtain the distribution of areas for orbits of all lengths one must 
integrate this against the exponential distribution of lengths discussed above. The 
leading behavior of this integral may be obtained by the method of steepest des- 
cent [85]: one finds the exponential form for P(|0|) with a c i ~ Vldd/Aa [89]. 
Numerical results shown in Fig. 10 clearly confirm this result for strongly chaotic 



scattering. 



Assuming fully chaotic scattering we insert the exponential form into Eq. (4. 1 9) 
and obtain a Lorentzian B-dependence, 

6R D {B)= * Y2 > < 4 - 2 °) 

1 + (2P/a c/ 0or 

where the magnetic field scale is set by the flux through the typical area a~ t l 
in units of the hc/2e (note the now-familiar appearance of the superconducting 
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Fig. 10. Distribution of classical effective areas in scattering from circle (top) showing power- 
law shape and strongly-deformed circle (bottom) showing exponential form characteristic of chaotic 
scattering. 

flux quantum in weak localization). So, as in the disordered case, for ballistic 
conductors with chaotic scattering, a weak magnetic field is sufficient to elimin- 
ate the enhanced backscattering. The effect this has on the conductance will be 
discussed in the next section. 

4.8. Problems with unitarity 

We have just shown using semiclassical theory that by application of a magnetic 
field B » a c /0o/2 the diagonal reflection coefficient is reduced from 271 to 1Z. 
By unitarity this means that Rod + T must increase by 11. However both Rod 
and the conductance (in units e 2 /h) T are of order N c and in order to calculate 
6T(B) we would need to evaluate them to order unity. But in section (4.3) we 
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Fig. 1 1. Energy-averaged weak localization magnetoconductance from the experiments of Keller ct 
al. obtained from data of the type shown in Fig. 9 compared to the prediction of Eq. (4.21) and of 
RMT. The parameter a used to scale the field axis was obtained from the conductance fluctuation 
power spectrum as discussed below, so there are no free parameters in the comparison. Inset shows 
one of the structures (referred to as the "stomach"). Also shown are numerical results from quantum 
simulations of the stomach. Adapted from [83]. 



were only able to evaluate T to order N c (T to order unity), so we are unable to 
give a semiclassical prediction for 6T{B) itself. We can make the following hand- 
waving argument: for an ideal chaotic system where there are no short trajectories 
all memory of the incident lead would be "lost" on average and all exit channels 
would be equally probable. Thus TZ = T = 1/2 and since Rod ~ R to order 
N' 1 , ST = 6 Rod = 1/4, so that 

1 + {2B/a cl (po) 2 

where we have now included for comparison with experiment the factor of two 
for spin. In fact exactly this result is obtained from a treatment using random- 
matrix theory [90]. However it must be admitted that a natural generalization of 
this argument to the disordered case gives gives a vanishing weak localization 
effect in quasi- ID [91], which is obviously incorrect. 

In Fig. 11 equation (4.21) is compared with experiment and excellent agree- 
ment is found. We remind you that there simply were no exact time-reversed pairs 
interfering in the semiclassical sum for T (and Rod) so this apparent failure of 
» the semiclassical method points to a deeper problem with the method in this con- 
text which is not yet resolved [72]. Clearly one aspect of this difficulty is that 
T and Rod must be evaluated to higher accuracy in the effective semiclassical 
parameter iV" 1 than Rd- 
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4.9. Ballistic conductance fluctuations 

We will not present here the full semiclassical theory used to describe the ballistic 
conductance fluctuations; it is described in detail in [72]. The basic technique is 
the same as the approach used above to treat WL. Here the objects of interest 
are the correlation functions of the conductance as a function of wavevector and 
magnetic field 



where 8T is defined by the deviation from the average value of section (4.3) 
in the absence of any symmetries, 8T = T — TkW/ir. The calculation of both 
correlation functions follows closely the work by Blumel and Smilansky [92] who 
computed semiclassically the correlation functions (in energy) of the S-matrix 
elements of chaotic systems. The correlation functions involve a sum over four 
set of classical paths 5, u, s', u' but with the terms s = u, s' — u' subtracted out 
since (as we have seen) they contribute to the average transmission. Wavevector 
averaging then picks out only the terms s = u' \u = s' which are only present in 
(8T a i,6T a b) (and not in terms of the form (8T a b8T c d)). In the resulting sum all 
phase factors are canceled except for those arising from the difference Afc, AB. It 
is natural then to compute the Fourier transform of these summations with respect 
to A/c, AB which we denote by C(L),C(S) respectively; these are just the power 
spectra of the fluctuations. Straightforward manipulations show [72] that they can 
be expressed in terms of the convolution of the distribution of path lengths with 
itself (in the case of C(L)) and of effective areas P(Q) with itself (in the case 
of C(S)). For the case of the length distribution (which is a simple exponential) 
its convolution is exponential with the same exponent, so C(L) oc exp[— *y c iL]- 
The area distribution is exponential in the absolute value of and so its self- 
convolution yields [70,72]: 



Note again, as in WL, a parameter a c i characterizing the average classical scat- 
tering dynamics determines a quantum property, in this case the power spectrum 
of the magnetoconductance fluctuations. A precise test of this semiclassical pre- 
diction is the following. First, an exact numerical solution for T{B) for a given 
chaotic cavity is obtained. The power spectrum of the magnetoconductance fluc- 
tuations is fit to the semiclassical form (4.25) and the parameter a is extracted. 
Then P(Q) is computed by Monte Carlo simulation of classical scattering in the 
same shape cavity and the parameter a c i is extracted. According to the semiclas- 
sical theory these two parameters should agree when N c » 1. 



C{Ak) = (8T(k + Ak)8T{k)) k 
C{AB) = {8T{k, B + AB)8T(k, B)) k 



(4.22) 
(4.23) 



C{S) a e 



a " |s| (l + a ci |S|). 



(4.24) 
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Fig. 12. Ratio of a obtained from quantum power spectra to that obtained from classical area distri- 
bution (see Fig. 10) for a range structures (see insets) vs. a c iA where A is the area of the structure. 

One finds [70 ,72] very good agreement over two orders of magnitude variation 
in a c i (see Fig. 12). Several experiments [64.65] have analyzed their magneto- 
transport data in exactly this manner and found fair but not impressive agreement 
with theory (in the experimental case the potential defining the cavity is imper- 
fectly known, and disorder [88] is apparently playing some role [891). 

4.10. Distinguishing regular from chaotic scattering 

If the scattering is not chaotic (if e.g. the cavity is circular or polygonal) then the 
distribution of effective areas P(0) is not exponential, but typically is a power 
law as shown in Fig. 10 [72,89] (because there are many different time scales 
for escape in an integrable system, depending on the invariant torus in the neigh- 
borhood of the scattering trajectory). In this case one still expects WL and CF 
but with statistical properties different from those derived above assuming hard 
chaos. For example the dependence of the WL magnetoconductance on magnetic 
field should deviate from the Lorentzian of Eq. (4.21) and the power spectrum 
of the fluctuations should deviate from Eq. (4.24). The observability of these 
deviations is perhaps the most interesting aspect of this theory because it implies 
that through purely quantum measurements one can determine the nature of the 
classical dynamics of a quantum system. Although the arguments leading to Eq. 
(4.2 1 ) must be modified for integrable systems, we can estimate the WL lineshape 
in such systems by simply taking P(O) ~ _/? in Eq. (4.19) ((3 « 3.3 for the 
circle). Since Eq. (4.19) is (up to scale factors) just the Fourier transform of P(0) 
with respect to B one expects a cusp-like singularity in 8T(B) as B —> 0, since 
the Fourier transform of a pure power law must have a singularity in the some 
derivative as the transform variable tends to zero. 




B (Gauss) 

Fig. 13. Weak localization magnetoresistance for (a) 48 stadium cavities, and (b) 48 circle cavities 
measured at T = 50 mK. The circular (integrable) cavities show the remarkable cusp-like shape sug- 
gested by the semiclassical theory, whereas the stadium (chaotic) cavities do not. Reprinted from [81]. 

In Fig. 13 is shown WL magnetoresistance data from an experiment by Chang 
et al. [81] on ensembles of ballistic microstructures (measured essentially in par- 
allel to average out the fluctuations). The ensembles differ only in the shapes of 
the electron cavities: stadium-like (top) and circular (bottom). Note the striking 
difference in the lineshape, the chaotic (stadium) ensemble showing smooth be- 
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havior as B — > as compared to the regular(circular) ensemble showing a sharper 
cusp-like behavior as expected from our rough argument. 

Well before this WL experiment Marcus et al. [64,79], motivated by similar 
reasoning, had studied the power spectra of magnetoconductance fluctuations in 
stadium and circular cavities. Since trapped trajectories in a circular cavity do not 
diffuse but instead circulate in a single sense one expects more enclosed flux on 
average and hence more high-frequency power over the same field range in com- 
parison to the stadium structure. This effect was observed and although it is not 
as dramatic in the data as the effect shown in Fig. 13, together these two exper- 
iments provide solid confirmation that quantum interference effects in ballistic 
microstructures actually can be used to distinguish quasi-regular from chaotic 
electron dynamics. It should be noted however that these experiments are diffi- 
cult and other experiments on integrable cavities with polygonal shape have failed 
to detect a significant deviation from the chaotic behavior [83.93J. 

4. 11. Other aspects of chaotic ballistic transport 

Another theoretical approach to chaotic ballistic transport is via random-matrix 
theory (RMT) [94]. This approach assumes that the junction is described by a 
hamiltonian from the Gaussian ensembles strongly coupled to the leads [90,95] 
or equivalently that the S-matrix of the junction is described by the Dyson circular 
ensemble [96.97]; it can calculate the amplitude and lineshape of the WL effect 
as well as the variance of the conductance fluctuations. This approach has no 
problem with unitarity to the required order and yields the value 8Twl = 1/4 
[95-97] as well as the Lorentzian lineshape [90], in agreement with Eq. (4.21) 
(but without the hand-waving to get the amplitude). It also predicts conductance 
fluctuations with a variance ((ST) 2 ) = 1/8/? where (3—1,2 depending on the 
presence or absence of time-reversal symmetry. It cannot address at all the non- 
chaotic junctions to which the RMT ensembles do not apply. 

It has now been understood [14,98,99] that although disordered microstructures 
are certainly classically chaotic their statistical properties differ from those pre- 
dicted by the standard Wigner-Dyson ensembles. This difference comes about 
essentially because a disordered system is described by a product of random 
matrices [100], which leads to a different form for the eigenvalue repulsion in 
the S-matrix [99]. These difference are relatively small [99] when comparing 
the ballistic WL and UCF effects to quasi-one-dimensional disordered systems 
(in the metallic regime), but become greater in two and three dimensions or when 
strong disorder leads to localization. Thus quantum interference effects in ballistic 
chaotic junctions are not identical to those of disordered systems. The ballistic 
junctions may be regarded as the "quasi-zero-dimensional" limit of a disordered 
quantum conductor. 
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Another interesting system is the ballistic quantum dot [62,63,66] in which 
the electron cavities we have been discussing are separated from the ieads by 
tunnel junctions. In this case the Coulomb charging energy [67] associated with 
the dot determines the spacing of conductance resonances, but we have proposed 
that the resonance amplitude is controlled by quantum interference effects [68]. 
The probability distribution for the amplitude may be calculated from standard 
RMT [68,101]. Recent experimental results find excellent agreement with this 
theory [102]. 

4.12. Summary 

Standard transport theory has to be extended to treat mesoscopic systems since 
the transport coefficients depend on sample and measurement geometry and on 
the particular impurity configuration in the case of disordered conductors. Quant- 
ities must be studied which are meaningful for a given mesoscopic system and 
not just for the average over an ensemble. The scattering approach of Landauer, 
Biittiker and Imry is a powerful approach to this end. It was shown above to 
be rigorously equivalent to quantum linear response theory applied to an infinite 
system of non-interacting fermions. In the case of ballistic junctions it naturally 
describes both classical magnetic and geometric effects and quantum interference 
effects via semiclassical theory of the transmission coefficients. Theory and re- 
cent experiments on these systems indicate that it is possible to distinguish chaotic 
from regular electron dynamics based solely on the statistical properties of the 
quantum interference phenomena in transport. 
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